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Abstract 

Anti-dynamo modes are usually found in spheromaks plasma devices experiments due to the 
fact that Cowlings anti-dynamo theorem is naturally applied to axisymmetric devices and flows. 
In this paper full consideration is given to the existence of slow dynamo modes in the case of 
compact Riemannian plasma devices without boundaries, such as tokamaks, stellarators and 
torsatrons. It is shown that a perturbed untwisted flow given by a decaying mode magnetic 
field is able to generate a slow twisted dynamo plasma flow where the unperturbed flow is also 
a steady flow. When the slow dynamo limit is obtained for high Reynolds magnetic numbers 
Rm, the unperturbed plasma flow achieve equilibrium as in spheromaks. A Riemann metric 
of a twisted thick magnetic flux tube with a very low aspect ratio is used in the computations. 
The data of aspect ratio and other data obtained from the Brazilian spherical tokamak at 
the National institute for space research (INPE) are used to obtain a numerical estimate 
for the maximum magnetic growth rate of the magnetic field of the slow dynamo action as 
7imax ~ ^^^'bT ^ toi'oidal initial magnetic field Bq ^ OAT and an aspect data tokamak 
of 1.5. For a tokamak internal radius of a ^ 0.2m one obtains a maximum growth rate for 
the slow dynamo as jimax ~ l-3x^, where vig is the poloidal perturbed diffusive flow. 
The safety factor is given by g ~ — 1 which is g < 1 implying stability of spherical tokamak 
plasma. Slow dynamos have been recently given in the literature as an example of Vishik's 
anti-fast dynamo theorem [Phys Plasmas 15 (2008)], which can also been tested in plasma 
experimental devices. Key-words:anti-dynamo theorems, slow dynamos. 
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I Introduction 



Recently Bellan [1] has been discussed the existence of anti-dynamo modes in plasmas sphero- 
maks. The stability has been given in terms of the safety factor conditions of instability {q > 1) 
or stability {q < 1). The Cowling's anti-dynamo theorem [2], which demands the existence of 
non-axisym metric magnetic flows for the dynamo existence, can be circumvented in sphero- 
maks as, even in axisymmetric spheromak devices, the existence of non-axisymmetric fluids 
and fields may exist, which per se, may guarantee dynamo action. Cowling [2], have show us 
how to recognise dynamo action in most types of magnetized plasma flows. The geometry of 
spheromaks, may be obtained by taking a very low aspect ratio, which is the ratio between 
the outer torus radius R and the cross-section a, or in the language of Riemannian twisted 
magnetic flux tube, the torus would be tick. Just as an example, in tokamaks in general the 
aspect ratio can be built in plasma laboratories, as A = ^ = 3, whereas in the Brazilian 
spherical tokamak [3], the ratio A — 1,5. Actually magnetic plasma stability is favour for 
aspect ratios of A <2. Indeed, the magnetic topology of the field lines in the spherical toka- 
maks are highly torsioned at its center, which allows us these fields to reach values as high as 
4T. These torsion central lines are used to concentrate a strong magnetic field. This leads us 
naturally to use the Riemannian geometry of twisted magnetic flux tube [4, 5] to model these 
spherical tokamaks, provided one considers only thick flux tubes where the aspect ratio is low. 
Another advantage of its use is that this leads easily to plasma stability. Another advantage is 
that according to Schuessler [6] astrophysical magnetic dynamos may be modelling by the flux 
tubes, as happens in the sun and other stars. Also another dynamo action mechanism called 
stretch-twist and fold (STF) as created by Vainshtein-Zeldovich [7] may naturally appear in the 
magnetic lines of the spherical tokamaks, since a fundamental ingredient of twist is the Frenet 
torsion of its lines in the dynamo flow [8]. From the mathematical side, using a dynamo model 
in Riemannian space has a long and solid tradition in toroidal maps as shown previously by 
Arnold et al [9] and in more recent papers by the author [5] to model plasma toroidal devices 
and astrophysical conformal dynamos. Another fundamental ingredient for dynamo existence, 
the folding, is a quantity which can be associated with Riemann curvature [10] to provide the 
doubling of the magnetic field intensity by a repetition process that may guarantee dynamo 
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action. In this paper untwisting flows are initially perturbed to yield slow dynamo action by 
using the spherical tokamak INPE data. In the presence of diffusion the relation between the 
growth rate of the perturbed magnetic field and the previous unperturbed one in terms of dif- 
fusion coefficient or resistivity. This allows us to show that a slow dynamo is obtained since the 
growth rate 7 vanishes when diffusion 77 0. Another Riemannian model for dynamo action 
has been recently obtained by Shukurov et al [11] to model small-scale dynamos such as the 
Perm toroidal experimental dynamo [12]. Another interesting distinction between the present 
model and the one addressed in reference (1112) is that in the case presently considered the 
perturbation of the magnetic field is not stationary. Recently a modification of Shukurov et al 
proposal was presented by the author can be considered as a thick Riemannian dynamo as well. 
Basic difference between our models and Moebius dynamo flow one proposed by Shukurov et 
al is that theirs is a numerical simulation and ours. Much earlier Mikhailovskii [13] has used 
a non-diagonal Riemann metric to describe tokamaks and investigate the plasma instabilities. 
The paper is organized as follows: In section II the absence of diffusion is shown to lead to 
non-dynamos or marginal dynamos. In section III diffusion is turn on and the perturbation 
is shown to lead to a slow dynamo spherical tokamak. Discussions and future prospects are 
presented in section IV. 
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II Riemannian plasma non-dynamo devices 



This section addresses the mathematical formalism of a general and thick twisted Riemannian 
magnetic flux tube and show that a marginal dynamo is obtained when a stationary unperturbed 
model is used when the magnetic field is given along the magnetic lines, which is B = 
B{s)t, and no slow dynamo is obtained. Here one considers that the growth rate of magnetic 
perturbed and unperturbed fields given respectively by Bi and Bq are independently given by 
7i and 70 where the magnetic fields are proportional to 

|Bo|ftie^°* (II. 1) 

and 

|Bi|«e^i* (II.2) 
where the perturbation process is given by 

B = Bo + Bi (II.3) 

where |Bo| >> |Bi|. Here t is part of the Frenet vector frame along the curve coordinate 
given by s-parameter. The complete Frenet frame is given by (t, n, b) where vectors n and 
b are the vectors that lay in the orthogonal plane to the vector t along the magnetic axis of 
the toroidal device. In this section one shall consider the diffusionless case, and show that this 
leads to a non-dynamo or at best to a marginal dynamo where 70 vanishes. But before digging 
into the physics of the problem let us take a moment to consider the Riemannian geometry 
of flux tubes as given for the first time in the context of solar plasma physics by Ricca [4]. 
The magnetic flux tube coordinates (r, ^r,s), is also used in plasma toroidal devices called 
tokamaks. Since folding processes in flux tubes can be represented by the Riemann curvature 
tensor, destructive folding that leads to non-dynamos in diffusive media, can be obtained by 
the vanishing of folding or vanishing of the Riemann curvature tensor. General flux tube 
Riemannian metric is 

dso' = dr'' + r'^de'' + K''{r, s)ds'' (II.4) 

The thin Riemann-flat in twisted magnetic flux tube metric is obtained by the constraining the 
relation K"^ :— {l—rK{s)cos9) to one. This is obtained as coordinate r approaches zero. This 
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happens in the neighbourhood of the torsioned flux tube axis. Coordinate 9{s) is one of the 
Riennannian curvilinear coordinates (r, 6r, s) and 6{s) — Or — J T{s)ds. The scalar function 
T represents the Frenet torsion. The thin tube metric is 

dso^ = dr^ + r^dOn^ + ds^ (II.5) 

The torsion term is responsible for the twist of the tube. Solar flux tubes are closed in the 
inner parts of the Sun , and then the tubes can be considered as compact Riemannian manifold 
without boundaries. Riemann gradient compact operator is given in general diffusive substrate 
by 

V = erdr + -de + t^d^ (II.6) 
r K 

while general self-induction equation is 

dtB = (B.V)v + 77AB (II.7) 

where in this section the resistivity 77 = in the ideal plasma case. Throughout this paper, 
the magnetic field is strictly confined along and inside the tube, which allows us to simplify 
the computations by considering that B^. — Q and that dgBg — 0. In the case consider here 

dtB = 78 - To^on (II.8) 

whose extra term is a non-inertial term similar to one that is introduced into a inertial frame 
by the use of curvilinear coordinates or Coriolis force in the frame. Therefore the diffusionless 
self-induction equation is given by 

[7ot - To'njSo = B^v^dstx (II.9) 

where due to the highly torsioned character of the internal spherical tokamak, one has used 
the helical hypothesis of circular helices where the torsion equals the Frenet curvature and 
are constants. By comparison of the both sides of equation (11.9) one obtains 

70 = (11.10) 

and 

= -To (11.11) 
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Here the Frenet equations 

dtt = nn (11.12) 

dtn = -Kt + rb (11.13) 

and 

dth = -rb (11.14) 

have been used to obtain the above results vq is the constant flow and the RHS of the equation 
(11.9) represents the stretching of the flow. Actually v = vot and the incompressibility condition 
of the flow 

V.v = (11.15) 

This implies actually that vq be constant. Thus from expression (11.10) one must conclude that 
no dynamo action is possible in for a constant modulus magnetic initial field in diffusionless 
media with a constant modulus stretching flow. Thus in the net section, we observe that 
by perturbing these magnetic fields with a non-stationary unsteady magnetic field a dynamo 
action is present but the dynamo is still slow. 

Ill Riemannian spherical tokamaks. slow dynamos 
and safety factor 

The above perturbation scheme described in the last section, is used in the above magnetic 
self-induced equation, where Bq and its perturbation Bi obey the following zero and first-order 
equations 

dtBo = (Bo.V)vo + 77AB0 (III.16) 

and 

dtB^ = (Bi.V)vo + (Bo.V)vi + 77AB1 (III.17) 
where A = is the Laplacian Riemannian operator given by 

A = [a^ + ia + ^ii±^a/i (111.18) 
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Applying the following perturbed magnetic field Bi into the above Riemannian Laplacian 
operator in the form 

Bi = Bieee + Bi,t (III. 19) 

where the dynamo mode m of the spatial toroidal coordinate-s appears as 

B^, = 6,6^1*+^"^^ (III.20) 

where m e Z, where Z is the field of integers numbers. Here bs is a constant with magnetic 
field dimensions. Note that the poloidal magnetic field perturbation Big is given by the same 
expression with the only difference that the constant bg of the toroidal field is now replaced by 
bi0. Note that the above choice seemingly hide the coordinate-^, but actually this is present 
in the ms exponent since the poloidal angular coordinate inside the tokamak, does depend on 
coordinate-s as ^ = 9{s) as contained above. The reason one is not using a radial mode in the 
magnetic field exponent is that one shall be considering here only compact spherical tokamak 
surfaces, where the radial dependence of the magnetic fields does not exist and the internal 
cross-section radius of the tokamak is a. Bellow one shall show that the slow dynamo mode 
m vanishes. Let us now display the expressions for the Laplacian of the initial field Bq = Bgt 
which by the solenoidality property 

V.Bo = (III.21) 

implies that Bq does not depend on coordinate-s, since however the tangent vector t = t(s), 
the vector field Bq is a non-uniform unsteady magnetic field. Due to the absence of radial 
dependence in the magnetic field over the tokamak Riemannian compact surface without 
boundaries the Laplacian of Bq simplifies to 

ABo=[fl4^5o(-t + n)] (III.22) 

From this expression, after some algebra, the complete diffusive self-induction equation results 
in the following growth rate 

(1 + cosef 

[acosO)'' 

Note that at this time, though the unperturbed magnetic field growth rate does not vanish 
the situation is even worse here, cause 70 is negative and this means that the initial toroidal 
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magnetic field decays as happens with some primordial magnetic fields in the universe [14]. 
The other equation yields the value of torsion r as 

t{s) « ^ (III.24) 

This result seems to be rather interesting since it shows that folding and twisting of the 
magnetic axis depends on a straightforward manner from the Let us now show that this does 
not happens with the perturbed field Bi growth rate 71. The solenoidal property of the 
perturbed field is 

V.Bi = (III.25) 
Application of the appropriate magnetic field into this equation yields 

^ = cose (III.26) 
Bis 

and 

im(cos^e - l)Bis = (III.27) 

this las equation leads to the constraint m = for the dynamo action mode if it exists at all. 
Substitution of this mode into the Riemannian Laplacian expression for the perturbed field 
reduces it to 

This yields the following growth rate as 

ris Bo . ^ 2rf{l + co?ff,sin2e , /xx. 

Note that in the limit 77 ^ the growth rate 71 ^ which characterizes the slow dynamo 
model. Thus one may say that the mode m = represents a slow dynamo mode. More general 
modes may be found which may represent a fast dynamo action in the spherical tokamak. Let 
us now consider the data involved in the INPE brazilian spherical tokamak experiment, in the 
case of maximum value of the growth rate ^imax, by taking into account the case when the 
magnetic Reynolds number Rm is high or when the dynamo flow is highly conductive. In this 
case the diffusion rj is small but finite and in this case the term of second order in -q may be 
dropped and finally the last expression reduces to 

llmax ~ ^ (III.30) 

3 a Bis 



Note that by using the INPE spherical tokamak data a = 0.2m and initial toroidal field 
Bq — OAT, this growth rate of slow dynanno is 

0.26 Vie , „ vw .TTT Q1^ 

llmax ~ S- ~ 1-3^ (III.31) 

a nis r>is 



here the torsion has been computed as 



To fti - fti 3.3m-^ (III.32) 
a 



Finally let us compute the safety factor q of the tokamak as 

d$ 1 ds 



de R rds 



(III.33) 



where we have used the expression above relating 9 and torsion integral / T{s)ds. Since Frenet 
curvature is given by 



These two last expressions yields 



«o = ^ (111.34) 



q = (III.35) 

To 



which in the circular helix case is given by g = — 1 < 1 and thus guarantees the tokamak 
plasma stability. 
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IV Conclusions 



Anti-dynamo modes in several plasma devices have been known in the literature. In this 
paper it is shown that slow dynamo modes can be obtained in non-turbulent flows diffusive 
plasma media. It is also shown that the absence of diffusion forbides the presence of a fast 
dynamo action in the case of initially toroidal flows where initial flows are aligned with the 
magnetic field. Data from the brazilian spherical tokamak operating in INPE is given to 
estimate the value of the growth rate of perturbed flows magnetic field. It is shown that 
this depend upon directly of the inverse of the magnetic Reynolds number which displays an 
explicitly slow dynamo behaviour. The Riemannian flux tube model in the thick case lead us 
to transform a toroidal tokamak into a spherical one, and the growth rate of the perturbed 
field is also proportional to the perturbed flow. When the perturbed flow vanishes the growth 
rate vanishes as well and no dynamo action is found whatsoever. Several other modes may be 
investigated in spherical tokamaks with the hope fast dynamo action may be found in future 
experiments. For example a more complicated model to INPE tokamak would be given by a 
tube in the center of the spherical tokamak where the torsion of the magnetic field could be 
confined around a Riemannian torus. 
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